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ABSTRACT 

We present arcsecond-resolution Submillimeter Array (SMA) polarimetric observations of the 
880 Ilia continuum emission from the protoplanetary disks around two nearby stars, HD 163296 and 
TW Hydrae. Although previous observations and theoretical work have suggested that a 2-3% polar- 
ization fraction should be common for the millimeter continuum emission from such disks, we detect 
no polarized continuum emission above a 3<j upper limit of 7mJy in each arcsecond-scale beam, or 
< 1% in integrated continuum emission. We compare the SMA upper limits with the predictions from 
the exploratory Cho & Lazarian (2007) model of polarized emission from T Tauri disks threaded by 
toroidal magnetic fields, and rule out their fiducial model at the ~ lOcr level. We explore some poten- 
tial causes for this discrepancy, focusing on model parameters that describe the shape, magnetic field 
alignment, and size distribution of grains in the disk. We also investigate related effects like the mag- 
netic field strength and geometry, scattering off of large grains, and the efficiency of grain alignment, 
including recent advances in grain alignment theory, which are not considered in the fiducial model. 
We discuss the impact each parameter would have on the data and determine that the suppression of 
polarized emission plausibly arises from rounding of large grains, reduced efficiency of grain alignment 
with the magnetic field, and/or some degree of magnetic field tangling (perhaps due to turbulence). 
A poloidal magnetic field geometry could also reduce the polarization signal, particularly for a face-on 
viewing geometry like the TW Hya disk. The data provided here offer the most stringent limits to 
date on the polarized millimeter- wavelength emission from disks around young stars. 
Subject headings: circumstellar matter — planetary systems: protoplanetary disks — polarization — 
stars: individual (HD 163296, TW Hydrae) 
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1. INTRODUCTION 

The magnetic properties of circumstellar disks are cen- 
tral to a wide range of physical processes relevant for 
planet formation. Dust and gas transport and mix- 
ing (e.g. Ciesla 2007), meteoritic composition (e.g. Boss 
2004), disk chemistry (e.g. Semenov et al. 2006), and 
the migration of planetary embryos through the disk 
(e.g. Chambers 2006) are all thought to be inffuenced 
by magnetohydrodynamic (MHD) turbulence. But per- 
haps the greatest impact of a magnetized disk is that 
MHD turbulence can provide the source of viscosity that 
drives disk evolution. Since the seminal work by Lynden- 
Bell & Pringle (1974), the photospheric excess and vari- 
ability exhibited by pre-main sequence stars have been 
attributed to an accretion disk. The viscous transport 
mechanism that supports the accretion process can also 
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explain many aspects of the time evolution of circumstel- 
lar disks (Hartmann et al. 1998), and by extension can 
help to constrain the physical conditions and timescales 
relevant for planet formation. However, there are re- 
markably few observational constraints on the magnitude 
and physical origin of viscosity in circumstellar disks. 

As conjectured by Shakura & Sunyaev (1973), tur- 
bulence can provide large enough viscosities to account 
for accretion and disk evolution on the appropriate 
timescales. The mechanism most commonly invoked as 
the source of this turbulence is the magnetorotational in- 
stability (MRI), in which magnetic interactions between 
fluid elements in the disk combine with an outwardly 
decreasing velocity field to produce torques that transfer 
angular momentum from the inner disk outwards (Balbus 
& Hawley 1991, 1998; see also Velikhov 1959 and Chan- 
drasekhar 1960). Indeed, it is unlikely that turbulence in 
an unmagnetized, azimuthally symmetric Keplerian disk 
can sufficiently redistribute angular momentum: mag- 
netic fields must be invoked to enable Shakura- Sunyaev 
viscosity (e.g. Balbus et al. 1996). The ionization frac- 
tion is likely high enough for magnetic coupling of ma- 
terial over much of the outer disk (see e.g. Sano et al. 
2000; Turner et al. 2007), and the observed Keplerian 
rotation of protoplanetary disks provides the requisite 
velocity shear. However, the magnetic field properties 
(strength and geometry) far from the central star remain 
unconstrained. 

Resolved observations of polarized submillimeter con- 
tinuum emission are uniquely suited to constrain the 
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magnetic field geometry - independent of disk structure 
- via the orientation of polarization vectors produced by 
dust grains aligned with the magnetic field (Aitken et al. 
2002). In the presence of an anisotropic radiation field, 
irregularly shaped grains with different cross sections to 
left and right circular polarizations of light can be spun 
up to high speeds by radiative torques (e.g. Dolginov 
1972; Dolginov & Mitrofanov 1976; Draine & Weingart- 
ner 1996)^. These spinning grains precess around mag- 
netic field lines, and ultimately align with their long axes 
perpendicular to the local magnetic field direction. Po- 
larized emission or absorption by these aligned grains can 
thus trace magnetic field structure in dusty interstellar 
media (see Lazarian 2007, and references therein). 

The first models of polarized emission from disks incor- 
porating the radiative torque alignment mechanism have 
recently been calculated by Cho & Lazarian (2007). Us- 
ing a two-layer Chiang et al. (2001) disk structure model 
threaded by a toroidal magnetic field (with circular field 
lines in the plane of the disk, centered on the star), they 
calculated the polarization emitted as a function of wave- 
length and position in the disk, incorporating emission 
and selective absorption mechanisms, but not scatter- 
ing. They predict a 2-3% polarization fraction at 850 /im, 
and note that grain alignment is particularly efficient in 
the low-density outer disk regions. At millimeter wave- 
lengths, dust grain opacities are low and optically thin 
thermal continuum emission primarily originates in the 
midplane where most of the mass is located. Polar imet- 
ric observations of millimeter- wavelength dust continuum 
emission therefore trace magnetic field geometry near the 
midplane in the outer disk, in regions where the magnetic 
field is strong enough for grains to become aligned and 
the density is low enough that grain spin-up is not im- 
peded by gas drag. 

The first attempt to observe polarized millimeter- 
wavelength emission from protoplanetary disks was made 
by Tamura et al. (1995, 1999). They used the James 
Clerk Maxwell Telescope (JCMT) to observe several 
young systems in the Taurus-Auriga molecular cloud 
complex - HL Tau, GG Tau, DC Tau and CM Aur - 
and reported tentative (~ 3a) detections of polarized 
millimeter- wavelength continuum emission from three of 
the four systems. The exception was GG Tau, for which 
they report a 2a upper limit of 3%. While the disks 
are unresolved in the 14" JCMT beam, the approximate 
alignment of the GM Aur and DG Tau polarization vec- 
tors with the known orientation of the disk minor axis 
is suggestive of a globally toroidal magnetic field struc- 
ture. However, finer resolution is required to confirm the 
magnetic field structure in the disks, and differentiate 
it from any potential contamination from an envelope 
or cloud material. DG Tau was followed up at a wave- 
length of 350 fim using the Caltech Submillimeter Obser- 
vatory by Krejny et al. (2009), and no polarization was 
detected with an upper limit of ^ 1%. They suggest that 
the decrease of polarization percentage relative to the 

^ More recent research in Lazarian & Hoang (2007) shows that 
in many cases rather than being spun up paramagnetic grains get 
slowed down by radiative torques, which means that the grains 
ahgned by radiative torques do not necessarily rotate suprather- 
mally. Nevertheless, the maximal rotational rate provides a useful 
parameterization of the effect of the radiative torques as discussed 
in detail in Hoang & Lazarian (2008). 



tentative 850 jam detection and the corresponding Cho 
& Lazarian (2007) prediction at 350 /im may be due to 
some combination of polarization self-suppression - effec- 
tively an absorption optical depth effect (see e.g. Hilde- 
brand et al. 2000) - or increased scattering at shorter 
wavelengths, which would produce a signal orthogonal 
to that expected for a toroidal magnetic field. In sum- 
mary, while the Cho & Lazarian (2007) predictions are 
consistent with the magnitude of the tentative JCMT 
detections, as discussed by Krejny et al. (2009), the pre- 
dicted polarization spectrum is inconsistent with the 350 
and 850 fim observations of DG Tau. 

In the absence of spatially resolved observations, the 
origin of the polarized emission in these systems remains 
unclear. While the position angle of the polarized emis- 
sion observed in the three systems in the Taurus- Auriga 
complex suggests association with the circumstellar disk, 
at least two of these sources (DG Tau and HL Tau) are 
fiat-spectrum sources host to jets and likely retain enve- 
lope material that could aid in generating a polarization 
signal (Kitamura et al. 1996; D'Alessio et al. 1997, e.g.). 
Tamura et al. (1995) suggest that the emission from HL 
Tau may arise from an interface region between the disk 
and a small envelope, and that the upper limit for GG 
Tau may be due to the lack of an envelope combined with 
weak, compact emission from the circumbinary ring. Ob- 
servations of the GM Aur system using the NICMOS in- 
strument on the Hubble Space Telescope indicate that it 
too may host a tenuous remnant outflow and envelope 
(Schneider et al. 2003). Nevertheless, with both theoret- 
ical predictions and observational evidence pointing to a 
2-3% polarization fraction at 850 /im in several T Tauri 
disks, resolved observations revealing the magnetic field 
geometry should be possible with current millimeter in- 
terferometers for bright disks. At ~ 1" resolution, such 
data would be well matched to the size scales at which 
the polarization fraction is expected to be the largest in 
the context of the Cho & Lazarian (2007) models. 

In order to test the Cho & Lazarian (2007) model pre- 
dictions and constrain magnetic field strengths and ge- 
ometries, we observed two nearby systems, HD 163296 
and TW Hya, with the Submilhmeter Array (SMA) 
polarimeter^. These targets were selected primarily 
for their large millimeter-wave fiuxes to maximize the 
expected polarization signal. Unlike the previously- 
observed Taurus targets, they are isolated from molecu- 
lar cloud material. HD 163296 has a total flux of 1.92 Jy 
at 850 /im (Mannings 1994), while TW Hya has a flux of 
1.45 Jy at 800 /im (Weintraub et al. 1989), predicting a 
total polarized fiux of ~40mJy in each system. Polar- 
ization of this magnitude should be observable with the 
SMA even if resolved across a few beams. HD 163296 
is a Herbig Ae star with a mass of 2.3 M© located at 
a distance of 122 pc (van den Ancker et al. 1998). It 
is surrounded by a fiared disk viewed at an intermedi- 
ate inclination of ~ 45°, observed to extend to at least 
500 AU in molecular gas and scattered light (Isella et al. 
2007; Grady et al. 2000), which has been extensively ob- 
served and modeled at millimeter wavelengths (Mannings 

The Submilhmeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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& Sargent 1997; Natta et al. 2004; Isella et al. 2007). 
TW Hya is a K star located at a distance of only 51 pc 
(Mamajek 2005; Hoff et al. 1998). It hosts a massive cir- 
cumstellar disk viewed nearly face-on at an inclination 
of 7° and extending to a radius of ~200 AU in molecular 
gas and scattered light (Qi et al. 2004; Roberge et al. 
2005). It is also a prototypical example of the class of 
disks with infrared deficits in their spectral energy distri- 
bution, known as "transition" disks. It has been shown 
to have a central deficit of dust emission extending out to 
4 AU (Calvet et al. 2002; Hughes et al. 2007), with a low 
mass accretion rate (Muzerolle et al. 2000) that may in- 
dicate clearing by a giant planet in formation (Alexander 
& Armitage 2007). 

We describe our observations of these systems in Sec- 
tion 2 and present the upper limits in Section 3. In Sec- 
tion 4.1 we describe the initial predictions generated by 
the Cho & Lazarian (2007) models and compare these 
predictions to the SMA observations. We then use these 
initial models as a starting point for an exploration of 
parameter space that seeks to describe how the different 
factors affect the predicted polarization properties (Sec- 
tion 4.2). We expand on these results by discussing the 
potential effects of physical mechanisms not included in 
the models (Section 4.3). In Section 5, we evaluate which 
physical conditions are most likely to contribute to the 
suppression of polarization relative to the fiducial model 
and summarize our results. 

2. OBSERVATIONS AND DATA REDUCTION 

Observations were conducted using the SMA polarime- 
ter, described in detail in Marrone & Rao (2008). The 
polarimeter uses a set of quarter-wave plates to con- 
vert the normally linear SMA feeds to circular polar- 
ization. By rotating the wave plate between two ori- 
entations separated by 90°, a single linear feed can 
be converted into either of the circular bases (left and 
right, or L and R). Since only one polarization (L or 
R) can be sampled on any antenna at any given time, 
full sampling of all four polarization states {LL, RR, 
RL, and LR) for all baselines (numbering N{N — l)/2, 
where N is the number of antennas) must be accom- 
plished by rotating the waveplates through a series of 
orientation patterns conveniently described by the two- 
state Walsh functions (see Marrone 2006). A full set 
of polarization states can be obtained for all baselines 
in ~5 minutes using a series of short (10-second) inte- 
grations with the waveplates rotated through patterns 
described by Walsh functions for the appropriate num- 
ber of antennas. These 5-minute intervals of data are 
combined into quasi-simultaneous Stokes parameters for 
each baseline (/ = {RR + LL)/2, Q = {RL + LR)/2, 
U = i{LR - RL)/2, V = {RR - LL)/2), where / mea- 
sures the total intensity, V measures the intensity of 
circular polarization, and ^JQ^ ^\J^ gives the total lin- 
early polarized intensity, which can also be used to cal- 
culate the local fractional (or percent) linear polarization 

^JQ^ + U^jl. Note that the aligned states [RR and LL) 
separate the total intensity from the circularly polarized 
intensity, while the crossed states (LR and RV) provide 
information about the linearly polarized intensity. Of 
course, this represents an idealization. Two relevant non- 
ideal effects are (1) if the R and L gains are not perfectly 



matched, then some of the bright Stokes / flux can leak 
into Stokes V when the difference is taken between LL 
and RR^ and (2) instrumental "leakage" of left circularly 
polarized light through a nominally right circularly po- 
larized waveplate (and vice versa) can transfer Stokes / 
to the linear states Q and JJ . 

Polarimetric SMA observations of the HD 163296 disk 
at 880 /im wavelength were carried out in the compact 
configuration on 29 May 2008, and in the extended con- 
figuration on 12 July 2008. The weather was excellent, 
with the 225 GHz opacity below 0.05 both nights, reach- 
ing as low as 0.03 on the night of 12 July. The phases 
were also extremely stable on both nights. The projected 
baseline lengths spanned a range of 9 to 260 kA, provid- 
ing a synthesized beam size of 1^1 x 0^.^89 for the combined 
data set, using natural weighting (see Table 1 for details). 
The quasar 3c454.3 was observed for 2.5 hours through 
more than 90° of parallactic angle during its transit, in 
order to calibrate the complex leakages of the quarter- 
wave plates. The quasar J 1733- 130 was used to calibrate 
the atmospheric and instrumental gain, and the quasar 
J 1924-292 was observed at 45-minute intervals through- 
out the night to test the quality of the phase transfer 
from J1733-130 as well as the calibration of the quarter- 
wave plate leakage. Uranus was used as the flux calibra- 
tor, yielding a flux for J 1733- 130 of 1.62 Jy on the night 
of May 29 and 2.01 Jy on the night of July 12. Uranus, 
Callisto, 3c273, and 3c279 were included as passband cal- 
ibrators. 

Observations of the disk around TW Hya were con- 
ducted in the subcompact and extended conflgurations of 
the SMA during the nights of 25 January and 15 Febru- 
ary 2009, respectively. Due to the far southern declina- 
tion of TW Hya in combination with the stringent ele- 
vation limits imposed to avoid antenna collisions in the 
subcompact conflguration, the source was only observ- 
able for three hours on the night of 25 January, while a 
full six hours of observations were obtained on 15 Febru- 
ary. The weather was again excellent, particularly for 
the extended configuration, during which the 225 GHz 
opacity remained stable between 0.03 and 0.04 for most 
of the night. The projected baseline lengths in the final 
data set varied from 6 to 250 kA, providing a synthesized 
beam size of l'.'2x0'.'9 in the final data set (see Table 1). 
The instrumental polarization was calibrated by observ- 
ing 3c273 over 90 degrees of parallactic angle for three 
hours across its transit. The quasar J1037-295 was used 
as the gain calibrator, and 3c279 was observed once per 
hour to test both the quality of the phase transfer and 
the instrumental polarization calibration. The primary 
fiux calibrator was Titan, yielding a fiux for J1037-295 
of 0.64 Jy on the night of 25 Jan and 0.53 Jy on the night 
of 15 Feb. 3c279, 3c273, and J1037-295 were included as 
passband calibrators. 

The double sideband receivers were tuned to a central 
frequency of 340.75 GHz (880 /im) for the HD 163296 
observations and 341.44 GHz (877 /im) for the TW Hya 
observations, with each 2 GHz- wide sideband centered 
±5 GHz from that value. The correlator was config- 
ured to observe the CO (3-2) transition (rest frequency 
345.796 GHz) with a velocity resolution of 0.70 km s~-^. 
The data were edited and calibrated using the MIR soft- 
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ware package^, while the standard tasks of Fourier trans- 
forming the visibihties, deconvolution with the CLEAN 
algorithm, and restoration were carried out using the 
MIRIAD software package. For a summary of the ob- 
servational parameters, including the 3a upper limits in 
Stokes Q and U for the individual tracks and the com- 
bined data sets, refer to Table 1. The test quasars for all 
tracks were point-like and unresolved. We detect polar- 
ized emission from the test quasars independently in each 
data set with a polarization fraction of between 8 and 
12% and a direction consistent between lower and up- 
per sidebands, as expected for linearly polarized emission 
from quasars at these wavelengths (see Marrone 2006). 

3. RESULTS 

Figure 1 shows the Stokes / (unpolarized) visibilities as 
a function of distance from the phase center in the (u^v) 
plane, corrected for the projection effects due to the disk 
inclination as in Lay et al. (1997) (for the mathemati- 
cal definition of the abscissa, see Section 3.3 of Hughes 
et al. 2008). This is effectively the Fourier transform of 
the radial brightness distribution of the disk. Both the 
HD 163296 and TW Hya disks are well resolved with 
high signal-to-noise ratios. 

We detect no polarized emission, in the CO (3-2) line 
or 880 /im continuum, from the HD 163296 or TW Hya 
disks. The rms values achieved in Stokes Q and U for 
the combined (compact+extended) continuum data are 
2.4mJybeam~^ and 2.3 mJy beam~^, respectively, yield- 
ing a 3(7 upper limit in both data sets of 7mJybeam~^. 
Given the integrated Stokes / fiuxes of 1.65 Jy and 
1.25 Jy for HD 163296 and TW Hya (see Table 1), the 
Cho & Lazarian (2007) result predictions of 2-3% polar- 
ization at these wavelengths imply ^30-50 mJy of polar- 
ized flux. Even if the spatial distribution of polarized 
flux in the source differs from that of the unpolarized 
emission, we should be able to detect it given that we re- 
cover most of the Stokes I flux. Figures 2 and 3 compare 
the data with the fiducial model predictions (described in 
Section 4.1 below). The upper right panel of each figure 
displays the amount and direction of observed polarized 
flux for each source, while the bottom row presents con- 
tour maps for each of the individual Stokes parameters. 
The emission in Stokes Q and U (linear polarization), 
as well as in Stokes V (circular polarization), is consis- 
tent with noise. As noted in Section 2, since Stokes V is 
calculated as the difference between the measured right 
and left {RR and LL) circular polarization, the difficulty 
of calibrating the gains precisely enough to remove the 
inffuence of the bright Stokes / emission raises the rms 
value in this Stokes parameter relative to Stokes Q and 
C/, which are calculated instead from the crossed (RL 
and LR) polarization states. 

We can rule out calibration errors as the reason for 
the lack of polarized emission for three reasons: (1) The 
point-like test quasars and the similarity of the visibility 
profiles in Figure 1 with previous observations of these 
sources (see e.g. Isella et al. 2007; Hughes et al. 2008) il- 
lustrate both the success of the atmospheric and instru- 
mental gain calibration and the high sensitivity of the 
data set. (2) The detection of polarized emission from 
the test quasars in each of the data sets, with direction 

^ See http://cfa-www.harvard.edu/~cqi/mircook.html. 



consistent between sidebands, demonstrates the success 
of the instrumental leakage calibration. Furthermore, (3) 
several of the nights were shared with other SMA polar- 
ization projects and our solutions for the instrumental 
leakage between Stokes parameters for the eight quarter- 
wave plates were effectively identical to those derived by 
other observers, who successfully detect polarization in 
their targets. 

It is worthwhile to compare the rms noise achieved here 
with the limiting precision of the current SMA polarime- 
ter. Errors in alignment of the quarter- wave plates intro- 
duce instrumental "leakage" between Stokes parameters, 
allowing some of the ffux from Stokes / to bleed into 
the linear Stokes parameters. The instrumental leakage 
correction is quite small (< 3%) and can to a large ex- 
tent be calibrated by observing a bright point source as 
it rotates through 90° of parallactic angle. Nevertheless, 
the uncertainty of this correction under typical observ- 
ing conditions is ^0.2%, although this can be reduced 
to < 0.1% with parallactic angle rotation, provided the 
source polarization does not vary with time (Marrone 
2006). Given the 2mJybeam~^ rms noise from our ob- 
servations compared with the peak Stokes / fiuxes of 740 
and 470mJy beam~^ (~0.3%), our constraints on the po- 
larized flux are approaching the limit of what is achiev- 
able with the SMA polarimeter. 

It is diflicult to directly compare the observations pre- 
sented here with the Cho & Lazarian (2007) model pre- 
dictions and the Tamura et al. (1999) JCMT result. The 
2-3% polarization factor reported by both sources refers 
to the integrated emission over the entire spatial extent 
of the disk. Since the SMA spatially resolves the emis- 
sion from the disk, the limit on the percent polariza- 
tion varies with position across the disk. The emission 
structure is predicted to be quite complicated (Cho & 
Lazarian 2007), with the percent polarization increas- 
ing as a function of distance from the star, so there is 
no straightforward way to quote a single value for the 
percent polarization that can be easily compared with 
the data. By tapering the SMA visibilities with a Gaus- 
sian whose FWHM is equal to the diameter of the disk 
as measured by a truncated power law model (Hughes 
et al. 2008), we can simulate an unresolved observation, 
similar to the JCMT result from Tamura et al. (1999). 
Using this method, we place a 3cr upper limit of 1% on 
the total polarized flux from both disks. However, such 
an extreme taper severely down- weights the visibilities on 
the longest baselines, which still have very high signal- 
to-noise ratios (see Figure 1). This effectively neglects 
the majority of the data: when all of the spatially re- 
solved data are taken into account, the limits are much 
more stringent, but they must be compared with the 
more complicated predictions from the spatially resolved 
model. Furthermore, decreasing the resolution may be 
additionally detrimental in the case of more face-on disks 
like TW Hya: if the magnetic field is perfectly toroidal, 
then the resulting radial polarization signal will cancel to 
zero in a large beam, no matter how strong the emission. 
To give a rough estimate, the ^ 40mJy of integrated 
polarized flux predicted for a 2-3% polarization fraction 
resolved into a few beams might predict a peak flux den- 
sity of 20mJybeam~^, which is about lOcr above the 
~ 2mJybeam~^ noise in the data. However, a detailed 
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TABLE 1 
Observational Parameters" 





HD 163296 


TW Hya 




Compact 


Extended 


C+E 


Compact 


Extended 


C+E 


Parameter 


29 May 2008 


12 July 2008 




25 Jan 2009 


15 Feb 2009 




340 GHz Continuum 


Beam Size (FWHM) 


2^2x1^/3 


0'.'9x0'.'7 


1^0x0'.' 9 


4^7x2^0 


l^lxO'.'S 


1^2x0'.' 8 


PA. 


50° 


-8° 


7° 


-1° 


6° 


6° 


RMS Noise (mjybeam"^) 














Stokes 


11 


7.8 


5.7 


35 


6.8 


5.5 


Stokes Q & [/ 


3.8 


2.7 


2.4 


6.3 


2.4 


2.3 


Peak Flux Density (mjybeam"^) 














Stokes / 


996 


639 


739 


990 


450 


474 


Stokes Q U (3cr upper limit) 


<11 


<8.1 


<7.2 


<19 


<7.2 


<6.9 


Integrated Flux*^ (Stokes /; Jy) 


1.65 


1.79 


1.64 


1.24 


1.33 


1.26 


CO(3-2) Line 


Beam Size (FWHM) 


2(^2x1^4 


0^9x0(^7 


lyixiyo 


5'.'0xiy2 


1^1x0'.' 7 


1^2x0'.' 7 


PA. 


50° 


-8° 


17° 


-1° 


9° 


9° 


RMS Noise (mjybeam"^) 


25 


19 


15 


55 


14 


13 


Peak Flux Density (mjybeam"^) 














Stokes / 


6500 


2650 


3730 


1520 


1800 


3090 


Stokes Q & C/^ (3o- upper limit) 


<75 


<57 


<45 


<170 


<42 


<39 


Integrated Flux® (Stokes /; Jykms"^) 


110 


56 


95 


47 


13 


27 



^All quoted values assume natural weighting. 

^The rms in Stokes I is limited by dynamic range rather than 
sensitivity. 

•^The integrated continuum flux is calculated using the MIRIAD 
task uvf it, assuming an elliptical Gaussian brightness profile. 

^The rms for the line is calculated using a channel width of 
0.7kms-^ 

^The integrated line fiux is calculated by integrating the zeroeth 
moment map inside the 3<t brightness contours. 

comparison with the spatially resolved model predictions 
for each disk can give a more robust result. 

The highest signal-to-noise ratio in an image is 
achieved using natural weighting, which assigns each vis- 
ibility a weight inversely proportional to its variance. In 
the case of observations with the SMA polarimeter, the 
bandwidth and integration time are the same for each 
integration, so the visibilities are primarily weighted by 
system temperature. For this reason, we use natural 
weighting to generate all images presented here. Using 
the upper limits from the naturally weighted images, it 
is possible to make comparisons with predictions of the 
spatially resolved emission generated from the models of 
Cho & Lazarian (2007). We pursue this avenue of inves- 
tigation in the following section. 

4. ANALYSIS AND DISCUSSION 

The constraints on polarized millimeter wavelength 
emission from the disks around TW Hya and HD 163296 
are inconsistent with previous observational (Tamura 
et al. 1999) and theoretical (Cho & Lazarian 2007) work 
that suggested that a polarization fraction of 2-3% should 
be common among protoplanetary disks. The strin- 
gent limit on the polarization fraction, when investigated 
within the context of the Cho & Lazarian (2007) model, 
can provide clues to the physical conditions within the 
disk that may be responsible for the suppression of po- 
larized emission relative to the fiducial model prediction. 
We therefore use the code described in Cho & Lazarian 
(2007) to generate models of the emission predicted for 
the TW Hya and HD 163296 disks, using available obser- 
vational constraints on the disk properties as inputs, and 



compare these predictions to the upper limits from the 
SMA observations (Section 4.1). We then identify pa- 
rameters that are not well constrained by existing obser- 
vations, and which have the greatest effect on the polar- 
ized emission rather than unpolarized Stokes / emission. 
We vary these parameters and investigate their effects 
on the predicted polarized submillimeter emission. We 
infer the range of values over which the predictions are 
consistent with the observations as well as the interac- 
tions between parameters in the context of the models 
(Section 4.2). Finally, we investigate other effects not 
implemented in these models that may contribute to the 
suppression of polarized disk emission, and estimate the 
magnitude of their contribution (Section 4.3). 

4.1. Initial Models 

The Cho & Lazarian (2007) predictions employ a two- 
layered Chiang et al. (2001) model of the density and 
temperature structure of a protoplanetary disk, includ- 
ing a surface layer with hot, small dust grains and an 
interior with cooler, larger grains. Within this model, 
the elongated dust grains are allowed to align via the 
radiative torque mechanism with a perfectly toroidal 
magnetic field threading the disk. The dust grains are 
assigned a size distribution described by a power law 
dN (X r~^s^^^"(ir where is the number of grains of size r, 
and (7grain IS initially taken to be 3.5 (Mathis et al. 1977). 
The grains are also assigned a degree of elongation given 
by the ratio of long-to-short axis cross sections, 
where C± and C\\ are the polarization cross sections for 
the electric field perpendicular and parallel to the grain 
symmetry axis, respectively. The grain size is defined as 
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Fig. 1. — Real (top) and imaginary (bottom) Stokes / continuum visibilities for HD 163296 (left) and TW Hya (right) as a function 
of distance from the disk center in the (u,v) plane, corrected for projection effects due to the inclination of the disk to our line of sight. 
Error bars show the standard error of the mean in each 7kA bin. See Lay et al. (1997) for details of the deprojection process. The inset 
in the upper right of each plot shows the CO (3-2) moment maps in Stokes / for the two disks. The colors indicate the first moment 
(intensity- weighted velocity), and the contours show the zeroeth moment (velocity-integrated intensity) in intervals of 3 Jykms"-*^. The 
solid line marks the position angle of the disk as determined by Isella et al. (2007) and Qi et al. (2004). The size and orientation of the 
synthesized beam is indicated at the lower left of each moment map. 



r, such that C± = (l + a)7rr^ and C\\ = (1 — a)7rr^, where 
a parameterizes the degree of elongation. The ratio of 
the major and minor axes of the grain are then given by 
a/b= y/{l-\-a)/{l - a). The grain shape is assumed to 
be oblate as in Cho & Lazarian (2007), consistent with 
observational evidence described in Hildebrand & Drago- 
van (1995). The initial 2-3% polarization estimates are 
based on the parameters for the "typical" T Tauri disk 
investigated in Chiang et al. (2001). 

In order to generate a model prediction that can be 
compared with the upper limits from the SMA observa- 
tions, we adjust these parameters to reflect the best avail- 
able information about the grain properties and density 
structures in the disks around HD 163296 and TW Hya. 
The initial model inputs, with references, are summa- 
rized in Table 2. We use temperature and surface den- 
sity power law indices and outer radii derived from previ- 
ous SMA 345 GHz continuum observations (Hughes et al. 
2008). The temperatures are calculated from the stellar 
temperature and gas and dust densities and opacities as 
in Chiang et al. (2001), while the surface density is ad- 
justed to best reproduce the observed 880 /im continuum 
flux. The temperatures and surface densities calculated 
here are consistent with previously determined values 
(e.g. Isella et al. 2007; Hughes et al. 2008) to within a fac- 
tor of two. Variations can be attributed to differences in 
the vertical temperature structure and dust grain opac- 
ities assumed in the models. While these disk structure 
models do not precisely reproduce the observed bright- 
ness profile, they represent a reasonable approximation 
within which the parameters determining the polariza- 



tion properties of interest can be investigated. 



TABLE 2 
Initial Model Parameters 





HD 163296 


TW Hya 


Parameter" 


Value 


Ref. 


Value 


Ref. 


n (K) 


9330 


1 


4000 


2 


(R©) 


2.1 


1 


1.0 


2 


M* (Mo) 


2.3 


1 


0.6 


2 


P 


0.8 


3 


1.0 


3 


Q-inner (AU) 


0.45 


4 


4.0 


5,6 


ao (AU) 


200 


3 


60 


3 


rmax,i (^m) 


10^ 


4 


10^ 


7 


i 


46° 


4 


7° 


8 


d (pc) 


122 


1 


51 


9,10 


Eo (g cm-2) 


130 




170 





References. — (1) van den Ancker et al. (1998); (2) Webb 
et al. (1999); (3) Hughes et al. (2008); (4) Isella et al. (2007); (5) 
Calvet et al. (2002); (6) Hughes et al. (2007); (7) Wilner et al. 
(2005); (8) Qi et al. (2004); (9) Mamajek (2005); (10) HofT et al. 
(1998) 

"^Symbols as in Chiang et al. (2001): T*, i?*, and M* are stel- 
lar temperature, radius, and mass, respectively; p and Eq describe 
the surface density profile E(i?) = Eo(i?/lAU)~^; ao is the outer 
disk radius; and rmax,! is the maximum dust grain size in the disk 
interior. Additionally, we define ainner (disk inner radius), i (incli- 
nation), and d (distance). All parameters not listed here are equal 
to the fiducial input parameters from Chiang et al. (2001). 

We use the model routines to generate 400x600 pixel 
sky-projected images (i.e. with 6- and 8-milliarcsecond 
pixels for TW Hya and HD 163296, respectively, signifi- 
cantly more finely spatially sampled than the data) giv- 
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Fig. 2. — Comparison between the Cho Lazarian (2007) model and the SMA 340 GHz observations of HD 163296. The top row shows 
the prediction for the model at full resolution (left), a simulated observation of the model with the SMA (center), and the 2008 SMA 
observations (right). The grayscale shows either the total flux (left) or the polarized flux (center, right), and the blue vectors indicate the 
percentage and direction of polarized flux at half-beam intervals. The center and bottom rows compare the model prediction (center) with 
the observed SMA data (bottom) in each of the four Stokes parameters (/, Q, U, V, from left to right). Contour levels are the same in 
both rows, either multiples of 10% of the peak flux (0.9 Jy/beam) in Stokes / or in increments of 2a for Q, U, and V, where a is the rms 
noise of 2.4mJy/beam. The size and orientation of the synthesized beam is indicated in the lower left of each panel. 



ing the total continuum flux, percent polarization, and 
orientation of polarized emission at each position across 
the disk. The full-resolution model is shown in the upper 
left panel of Figures 2 and 3, although the lines indicating 
orientation have been vector- averaged in bins of several 
pixels for clarity of display. We then use the MIRIAD 
task uvmodel to sample the image with the same spatial 
frequencies as the SMA data. We invert the visibilities 
and image with natural weighting to create a simulated 
SMA observation of the disk model, shown in the top cen- 
ter panel of Figures 2 and 3. We also create simulated im- 
ages in each of the four Stokes parameters (center row), 
since the Stokes parameter images are most directly com- 
parable to the upper limits set by the observations. The 
model images show the distinctive quadrupolar pattern 
in Stokes Q and U predicted by the model for a toroidal 
magnetic field geometry, due to the radial orientation 
of the polarization vectors. The intermediate inclina- 
tion of HD 163296 creates an hourglass-shaped bright 
region along the disk minor axis, where the synthesized 
beam picks up emission from the highly polarized regions 
along the front and back of the outer disk, concentrated 
towards the disk center by the viewing geometry. This 
predicted morphology echoes the alignment of polariza- 
tion vectors with the minor axes of the disk observed by 
Tamura et al. (1999). With predicted peak Stokes Q and 



U fiuxes of 23 and 16 mJybeam ^, these initial models 
of polarized emission are ruled out at the lOa and To- 
level for HD 163296 and TW Hya, respectively, by the 
SMA upper limits. 

4.2. Parameter Exploration 

With the fiducial model prediction ruled out at high 
confidence, we turn to an exploration of the input param- 
eter space to provide information about the conditions in 
the disk that might be responsible for the suppression of 
polarized emission. We first identify several parameters 
that most strongly affect the polarization properties of 
the disk, without significant impact on the Stokes / emis- 
sion. In the Cho & Lazarian (2007) model, the radiative 
torque mechanism that spins up elongated dust grains 
along magnetic field lines is impeded primarily by gas 
drag in regions of high density. Since we normalize the 
surface density to reproduce the 880 /im fiux (for the as- 
sumed opacities and derived temperatures), we cannot 
vary this quantity. However, the degree of elongation of 
the dust grains, the threshold set within the model for 
grain alignment, and the dust grain size distribution are 
all important factors that affect the polarization proper- 
ties of the disk rather than the Stokes / emission. These 
parameters are discussed in greater detail in the following 
sections. 
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Fig. 3. — Comparison between the Cho &; Lazarian (2007) model and the SMA 340 GHz observations of TW Hya. The top row shows the 
prediction for the model at full resolution (left), a simulated observation of the model with the SMA (center), and the SMA observations 
(right). The center and bottom rows compare the model prediction (center) with the observed SMA data (bottom) in each of the four 
Stokes parameters (/, Q, U, V, from left to right). Contour levels are the same in both rows, either multiples of 10% of the peak flux 
(47mJy/beam) in Stokes / or at 2cr intervals for Q, U, and V, where a is the rms noise of 2.3mJy/beam. Symbols as in Figure 2. 



4.2.1. Grain Elongation 

The elongation of the dust grains is important both 
for the radiative torque and because the differing cross- 
sections parahel and perpendicular to the magnetic field 
allow the grain to emit polarized continuum emission at 
millimeter wavelengths. The fiducial model assumes a 
long-to-short axis cross-section ratio C±/C\\ =2.1, cor- 
responding to an axial ratio of 1.5:1 for oblate dust grains 
(for the relationship between cross section and axial ra- 
tios for different grain geometries, see e.g. Padoan et al. 
2001). Varying this ratio determines the radial extent 
of the disk over which the dust grains are aligned with 
the magnetic field, as well as how much polarized light 
is emitted from the disk: it effectively changes the effi- 
ciency of grain alignment and the emission cross-section 
of the grains. 

In order to obtain a quantitative description of the 
effect of grain elongation on the predicted intensity of 
polarized emission from the disk, we generate a series of 
models with different cross section ratios as described in 
Cho & Lazarian (2007) with initial parameters listed in 
Table 2. We then sample the model images with the SMA 
spatial frequencies, as described in Section 4.1 above, and 
compare the peak flux in Stokes Q and U with the 3cr 
upper limit from the SMA observations. Figure 4 plots 
the peak flux in the Stokes Q and U model images as a 
function of the dust grain cross section ratio. For com- 



parison, the shaded area marks the region of parameter 
space consistent with the 3<j upper limits from the SMA 
observations. The series of panels across the top of the 
plot show the model images in Stokes Q and [/, sampled 
with the SMA spatial frequencies, for three representa- 
tive values of the dust grain cross section ratio. From 
these maps, it is clear that the dust grain elongation acts 
primarily as a scaling factor for observations at this res- 
olution: the emission morphology does not change, but 
simply becomes stronger or weaker as the dust grains be- 
come more or less elongated. From the HD 163296 plot 
on the left and the TW Hya plot on the right, we can 
see that if the dust grain elongation were the only factor 
suppressing polarized emission from the disk, the grains 
would have to be quite round, with C±/C\\ < 1.2 — 1.3. 

4.2.2. Grain Alignment Griterion 

Another model input that is important for the polar- 
ization properties of the disk is the value at which the 
threshold for grain alignment via the radiative torque 
is set. In order to determine whether or not the dust 
grains are aligned with the magnetic fleld in a particular 
region of the disk, a comparison is made between the ro- 
tational kinetic energy imparted by the radiative torque 
and that imparted by random collisions with gas parti- 
cles in the disk. A useful parameterization is {oorad/^th)'^ , 
where corad and coth are the angular velocities of the grains 
due to radiative torques and thermal collisions, respec- 
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Fig. 4. — Peak continuum flux in Stokes Q and C/ as a function of dust grain cross section ratio for HD 163296 (left) and TW Hya 
(right). The top row shows the resolved emission in Stokes Q and U predicted for three values of the dust grain cross section ratio, sampled 
at the same spatial frequencies as the data. The grayscale indicates the intensity of emission relative to the peak flux of the data when 
the grain cross section ratio equals three, with white indicating positive emission and black indicating negative emission. Contours are 
[2,4,6,...] times the rms noise (2.4 mJy for HD 163296 and 2.3 mJy for TW Hya) with positive contours in black and negative contours in 
white. The plots below give the peak flux in the synthesized beam predicted by the models as a function of the grain cross section ratio. 
Stokes Q is plotted as a solid line while Stokes C/ is a dotted line. The three-sigma upper limit on the peak flux from the SMA observations 
is indicated by the gray region of the plot. The y-axis along the upper edge of the plot gives the dust grain axial ratio. All images and 
peak flux values assume natural weighting to minimize noise. 



tively. The radiative torques act to align grains with 
the magnetic field, while gas drag inhibits alignment and 
causes grains to point in random directions: the ratio 
(^rad/^th)^ therefore serves as a measurement of the ef- 
fectiveness of the radiative torque in aligning the grains 
with the magnetic field. This ratio will generally be high- 
est, and the grains most aligned, in the outer disk where 
the gas density is low. We therefore expect grains to be 
aligned in the outer disk, and oriented randomly in the 
inner disk. Since the value of (cc^rad/^th)^ varies with ra- 
dial distance from the star, the chosen threshold value for 
alignment effectively varies the radius at which grains be- 
come aligned with the disk magnetic field. The threshold 
is initially set so that grains are assumed to be aligned in 
regions of the disk where the kinetic energy imparted by 
the radiative torque is 10^ times greater than that im- 
parted by thermal collisions. We vary this threshold in 
order to study its effects on the polarization properties 
of the disk. 

Figure 5 shows the peak fiux predicted for Stokes Q 
and /7 as a function of the grain alignment threshold 
(^rad/^th)^, Compared with the 3<j upper limit from the 
SMA observations for HD 163296 (left) and TW Hya 
(right). It is clear that for both disks, the threshold 
would have to be set many orders of magnitude higher 
than the conservative initial value in order for the align- 
ment to be weak enough to account for the lack of a 
polarization signal. Indeed, in order for this to be the pri- 
mary mechanism suppressing the disk polarization, the 
threshold would need to be raised until alignment is per- 
mitted to occur only when the rotational kinetic energy 
imparted by the radiative torque is at least 5-7 orders of 
magnitude greater than that of gas grain collisions. This 



is most likely an unrealistically stringent constraint. 

It should be noted here that the approach to alignment 
in Cho & Lazarian (2007) requires revisions to account 
for recent advances in the quantitative theory of grain 
alignment. First of all, in the calculations of the ra- 
tio {oOrad/^th)i the simplifying assumption is made that 
the radiation seen by each grain is coming from a point 
source. In fact, the bulk of the radiation field originates 
as reprocessed starlight from neighboring regions of the 
disk, so although there should be an overall radial gradi- 
ent, it is better approximated by multipoles rather than a 
purely unidirectional signal. When the effects of this ra- 
diation structure are accounted for, the ratio (oorad/^th) 
can decrease by up to a factor of 10 (Hoang & Lazar- 
ian 2009, Figure 17). An additional decrease by another 
factor of ~10 may come from the fact that the overall 
direction of anisotropy is perpendicular to the assumed 
toroidal magnetic field in the disk (Hoang & Lazarian 
2009, Figure 17). This effect may be mitigated some- 
what in a clumpy disk, where local anisotropics will not 
necessarily be radially oriented and may even be aligned 
with the magnetic field. Taking both effects into account 
and squaring the ratio demonstrates that the kinetic en- 
ergy of the grains in their maximal state of rotation may 
be up to 4 orders of magnitude less than is assumed using 
ad hoc assumptions in the spirit of the old understand- 
ing of radiative torque alignment. An additional decrease 
comes from the fact that an appreciable portion of grains 
may be aligned in the so-called "zero-J" alignment point 
(Lazarian & Hoang 2007). Grains in this point are not 
perfectly aligned as assumed in Cho & Lazarian (2007), 
but instead will wobble, reducing the degree of alignment 
to only ~20% (see Hoang & Lazarian 2008). In addition. 
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while interstehar grains are always aligned with long axes 
perpendicular to magnetic field, larger grains in circum- 
stellar disks may not have efficient internal relaxation 
and can be occasionally aligned with long axes paral- 
lel to magnetic field (Hoang & Lazarian 2009). These 
factors can significantly decrease the observed degree of 
polarization expected from the circumstellar disks com- 
pared to the Cho & Lazarian (2007) estimate, making 
the predictions roughly comparable (to within an order 
of magnitude or so) to the SMA upper limits. 

4.2.3. Grain Size Distribution 

Cho & Lazarian (2007) emphasize the importance of 
the grain size distribution in determining the observed 
polarization properties of circumstellar disks. We fix the 
minimum grain size at rmin = 0.01 /im as in Cho & Lazar- 
ian (2007) and Chiang et al. (2001). Although growth to 
larger sizes may have occurred, the minimum grain size 
affects the millimeter- wavelength polarization properties 
in the context of the model only through the normaliza- 
tion of the total mass: increasing the minimum grain size 
to 1 /im (required to reproduce the 10 /im silicate feature 
from the inner disk; see e.g. Calvet et al. 2002) changes 
the predicted polarization by less than 0.1%, since it does 
not bring the density above the threshold value necessary 
to suppress grain alignment in the outer disk. Two as- 
pects of the grain size distribution that can be varied 
in the context of the model are the maximum grain size 
Tmax and the power law index ggrain, where the grain size 
number density goes as dN oc r~^srain^^ 

Observational evidence points to grain growth up to at 
least 1 mm in the HD 163296 disk (Isella et al. 2007) and 
1cm in the TW Hya disk (Wilner et al. 2005), without 
ruling out the possibility that grains have grown to even 
larger sizes (perhaps even planetary dimensions in the 
case of TW Hya; see Calvet et al. 2002; Hughes et al. 
2007). Since the surface density is chosen to maintain 
consistency with the observed 880 /im ffux in Stokes /, 
the number density of particles with sizes near 880 /im, 
which dominate the 880 /im flux, remains roughly con- 
stant regardless of the maximum grain size in the dis- 
tribution. Thus the effect of raising the maximum grain 
size in the distribution is primarily to introduce "invisi- 
ble" grains at sizes larger than 1 mm or 1 cm, which has 
no effect on the observable polarization properties (cf. 
Figure 7 in Cho & Lazarian 2007). However, adding 
mass at the large-grain end of the size distribution while 
keeping constant the mass in small grains has the effect 
of raising the total surface density of the disk. This is 
unrealistic for all but a small increase in maximum grain 
size, as the disk quickly becomes Toomre unstable and 
gravitational collapse or deviations from Keplerian ro- 
tation should rapidly become observable. While this is 
most likely an artifact of the assumed grain size distri- 
bution, it suggests that within the context of the model, 
grain growth is unlikely to be the mechanism suppressing 
the emission of polarized radiation. 

The power law index Qgrain 

controUing the relative pop- 
ulation of large and small grains in the disk is some- 
what more promising. In general, the polarized emission 
observed at a particular wavelength will tend to origi- 
nate primarily from dust grains smaller than the wave- 
length, while the unpolarized emission will be dominated 
by grains of roughly the same size as the wavelength. Be- 



cause dust grains of size ~880 /im are within the geomet- 
ric optics regime (27rr/A > 1, where A is the wavelength 
of observation, 880 /im), they do not contribute to the 
polarized emission predicted by the models. Most of the 
Stokes Q and U emission at these wavelengths originates 
from dust grains with sizes less than A/27r ~ 100 /im 
(Cho & Lazarian 2007), while most of the Stokes / emis- 
sion originates from grains with sizes similar to the wave- 
length of observation. The relative number of 100 and 
880 /im grains in the disk, determined by ggrain, there- 
fore plays a role in determining the amount of polar- 
ized emission observed. However, since the differences in 
grain sizes is not large, the power law index must change 
substantially before the effect on the polarization prop- 
erties becomes appreciable. Varying ggrain from 3.5 to 2 
changes the peak linearly polarized flux in the model by 
only 20%. Therefore, when comparing the SMA limits 
with the model predictions, the dust grain size distribu- 
tion has relatively little impact on the predicted polar- 
ization properties of the disks. 

4.2.4. Interactions Between the Parameters 

The analysis so far has explored individual model pa- 
rameters as though they were fully independent, deter- 
mining the range of values permitted by the SMA upper 
limit for each parameter separately. However, it is useful 
to understand how the parameters relate to one another 
in determining the polarization properties of the disk. 
Here we investigate relationships between pairs of the 
parameters considered above. 

We first study the relationship between dust grain elon- 
gation and the grain alignment threshold. As discussed 
in §2.3 of Cho & Lazarian (2007), the rotation rate of 
dust grains due to the radiative torque is a function of 
the peak wavelength of the radiation field and the dust 
grain size, with no explicit dependence on grain axial ra- 
tio. As described in Dolginov & Mitrofanov (1976), spin- 
up by the radiative torque mechanism is caused by the 
irregular shape of the grain, which gives it differing cross 
sections to left and right circular polarization; elongation 
does not necessarily favor either polarization basis. This 
is reffected in the table of timescales relevant for grain 
alignment in Lazarian (2007): neither the radiative pre- 
cession time nor the gas damping time depends on the 
grain axial ratio. The primary effect of the grain elon- 
gation in alignment is to decrease the Larmor precession 
time, which causes the spinning grains to align their ma- 
jor axes more quickly with the magnetic field lines (or, al- 
ternatively, decreases the critical magnetic field strength 
in a given region of the disk; see Section 4.3 below). We 
therefore do not expect much, if any, dependence between 
these variables. In order to test this expectation, we vary 
the dust grain cross section ratio and the grain alignment 
threshold for the HD 163296 disk. The model prediction 
of peak flux in Stokes U (which provides the most strin- 
gent limits when compared to the SMA data) are shown 
in Figure 6. The shaded gray region of the plot repre- 
sents the parameter space within which the model pre- 
diction is less than the 3a upper limit given by the SMA 
data, i.e., combinations of parameters consistent with the 
observational results. The contours show the predicted 
peak flux of the model in Stokes U for each combination 
of parameters: model predictions with greater polarized 
intensity are more strongly inconsistent with the obser- 
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Fig. 5. — Peak continuum flux in Stokes Q and t/ as a function of the threshold for grain ahgnment (see Section 4.2 in the text) for 
HD 163296 (left) and TW Hya (right). The top row shows the resolved emission in Stokes Q and U predicted for three values of the 
alignment threshold, sampled at the same spatial frequencies as the data. The grayscale indicates the intensity of emission relative to 
the peak flux of the data when the alignment threshold equals 10^, with white indicating positive emission and l3lack indicating negative 
emission. Contours are [2,4,6,...] times the rms noise (2.4 mJy for HD 163296 and 2.3 mJy for TW Hya) with positive contours in black 
and negative contours in white. The plots below give the peak flux in the synthesized beam predicted by the models as a function of the 
alignment threshold. Stokes Q is plotted as a solid line while Stokes C/ is a dotted line. The three-sigma upper limit on the peak flux from 
the SMA observations is indicated by the gray region of the plot. The y-axis along the upper edge of the plot gives the dust grain axial 
ratio. All images and peak flux values assume natural weighting to minimize noise. 



vational limits. Because of the assumption in the models 
that grains meeting the alignment criterion will become 
aligned with 100% efficiency, grain alignment and elon- 
gation are evidently only weakly coupled. 

Another potentially important relationship is that be- 
tween the grain size distribution and grain elongation. 
Little is known about the relationship between these vari- 
ables, since both are notoriously difficult to constrain ob- 
servationally. Nevertheless, if grains grow simply by ac- 
cumulating material evenly over their surface then they 
may naturally become more spherical as they become 
larger. Spherical grains emit less strongly polarized radi- 
ation than more elongated grains, so it might be expected 
that grain growth can suppress the emission of polarized 
light, even in cases where the alignment mechanisms are 
quite efficient (as expected for large grains, e.g. Cho & 
Lazarian 2005). Indeed, a corresponding inverse rela- 
tionship between grain size and polarization fraction has 
been observed in molecular clouds (e.g. Vrba et al. 1993). 
Given the observed growth to millimeter and even cen- 
timeter sizes within the disks around HD 163296 and 
TW Hya (Isella et al. 2007; Wilner et al. 2005), and the 
large (~100/im) sizes of the grains responsible for emit- 
ting most of the polarized radiation (see Section 4.2.3), it 
is perhaps plausible that the grains in these disks should 
have cross section ratios consistent with the values of 
1.2-1.3 constrained in Section 4.2.1 above. We know that 
this cannot be true everywhere in the interstellar medium 
(ISM): polarization at 850/im is observed in star-forming 
regions at much earlier stages (e.g. Girart et al. 2006), 
and far-infrared polarimetry indicates that grains with 
axial ratios a/b between 1.1-3 are common at sizes of tens 
of microns in molecular clouds (Hildebrand & Dragovan 



1995). However, a tendency towards spherical grains in 
T Tauri disks, even just at the low end of the distribu- 
tion inferred by Hildebrand & Dragovan (1995), should 
be able to suppress the emission of polarized radiation 
from the disk enough to bring the models within range 
of the observational constraints. 

We can test the plausibility of this degree of elonga- 
tion by modifying the discussion of grain growth based 
on turbulent coagulation in Vrba et al. (1993). If we 
assume that the grains in T Tauri disks originate exclu- 
sively from small, highly elongated grains in the ISM, 
e.g. with initial major axis = 0.1 /im and axial ra- 
tio ai/bi = 2 (Aannestad & Purcell 1973; Hildebrand 
& Dragovan 1995), then we can estimate how the ax- 
ial ratio changes with grain size. Neglecting asymmetric 
effects like collisional destruction, grain size might be 
expected to grow roughly evenly in all directions with 
the number of grain-grain collisions, A^, in such a way 
that the final grain size is simply aj = aiN'^^^. The 
change in any dimension of the grain, 5, is then given 
by (5 = aiN'^^^ — yielding a final minor axis size 
of 6/ = 6i + J, or a f/bf = af /{af — ai -\- bi). If ai 
= 0.1 /im and a/ = 100 /im, then af/bf = 1.001, sig- 
nificantly more round than the upper limit set by the 
SMA data. The timescale needed for grain growth to 
these (up to meter) sizes is of order 10^ years at a dis- 
tance of 50 AU from the central star (see e.g. Weiden- 
schihing 1988; Dullemond & Dominik 2005). This cal- 
culation is highly simplified and neglects complications 
like the evolution of conditions within the disk, shap- 
ing by grain-grain collisions (e.g. Dullemond & Dominik 
2005), and the complexity of the grain size distribution. 
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Yet the extremely spherical grains produced on relatively 
short timescales in this oversimplified scenario represent 
a lower limit to the grain elongation that perhaps sug- 
gests a scenario by which grains might have grown into 
shapes that are nearly spherical enough (with axial ra- 
tios of 1.2 rather than 1.001) to plausibly account for the 
suppression of polarized emission. 




KO K5 2.0 2.5 3.0 
Grain cross section ratio Cj^/Ch 



Fig. 6. — Detect ability of Stokes U continuum emission from 
the HD 163296 disk as a function of the dust grain cross section 
ratio (Section 4.2.1) and the threshold for dust grain alignment 
(Section 4.2.2). The gray regions of the plot represent portions 
of the parameter space that would be undetectable given the 3a 
upper limit from the SMA observations, while contours show the 
peak flux of the model for each set of parameters, beginning at the 
3(7 level (7.2 mJy) and increasing by intervals of 2a (4.8 mjy). The 
two parameters are only weakly degenerate. 

4.3. Other Effects 

In the previous section, we investigated the effects 
of those parameters considered in the Cho & Lazarian 
(2007) model. However, there are additional effects that 
may also play a role in suppressing polarized emission 
from the disk relative to the fiducial 2-3% prediction. 
Among these are the magnetic field strength, the geo- 
metric regularity of the magnetic field, and polarization 
due to scattering. 

4.3.1. Magnetic Field Strength 

The magnetic field strength plays a role in determining 
whether or not grains can become aligned via the radia- 
tive torque. If the magnetic field is above some criti- 
cal strength, grains will become aligned provided that 
the radiative torque can generate more rotational kinetic 
energy than thermal collisions. At low magnetic field 
strengths, grains are not expected to align with the mag- 
netic field at all. The critical magnetic field strength for 
alignment may be estimated by comparing the Larmor 
precession time tL with the gas damping time tgas- Fol- 
lowing Lazarian (2007) and using fiducial values for the 
magnetic susceptibility and dust grain density, alignment 
is possible when tL < tgas, or: 

B> 4.1x10-5^^^^^^^ (1) 



where B is the magnetic field strength in units of /iG, 
r is the grain size in cm, n is the gas density in units 
of cm~^, Td is the dust temperature in K, Tg is the gas 
temperature in units of K, and s is the ratio of minor 
to major dust grain axes. Using the power law mod- 
els of density and temperature derived in Hughes et al. 
(2008), it is possible to estimate these quantities for the 
regions of the outer disk probed by the SMA data. Tak- 
ing the values at disk radii equivalent to the spatial reso- 
lution of the data (-I'.'O, or 50 and 120 AU for TW Hya 
and HD 163296, respectively), and assuming equivalent 
gas and dust temperatures, we derive densities of sev- 
eral times 10^ cm~^ and temperatures of ~ 40 — 50 K. 
For the 10-100 /im grains contributing most of the polar- 
ized emission in the models, the critical magnetic field 
strength is of order 10-100 mG. 

This strength matches reasonably well with theoreti- 
cal expectations. Shu et al. (2007) developed a model 
of steady-state magnetized accretion disks that predict 
magnetic field strengths of order 10-100 mG on the spa- 
tial scales probed by the data. Wardle (2007) pointed 
out that Zeeman splitting of OH in molecular cloud cores 
and masers in star forming regions place a lower limit of 
~10mG on the magnetic field strength, which will likely 
be amplified by compression and shear during the pro- 
cess of collapse that forms the central star and disk. It 
should also be noted that the value quoted above should 
be taken as a lower limit, since superparamagnetic inclu- 
sions would significantly decrease the required magnetic 
field strength for alignment (Lazarian & Hoang 2008). 
The critical magnetic field strength required to align 
grains within the conditions of the model is therefore 
reasonable compared to theoretical expectations. We do 
not expect that the lack of polarized emission is due to 
extremely low magnetic field strengths. 

4.3.2. Geometric Regularity of the Magnetic Field 

The assumption that the field is toroidal arises from 
the supposition that the rotational motion of the disk 
has affected the magnetic field geometry. Yet for this to 
occur, the ionization fraction must be large enough that 
disk material and magnetic fields can interact. How- 
ever, this also implies that turbulent motions within 
the disk (perhaps even of magnetic origin) may tan- 
gle the magnetic fields locally, adding a random com- 
ponent to the ordered toroidal magnetic field. It is ex- 
tremely difficult to estimate the magnitude of such an 
effect without knowing both the ionization fraction and 
the magnitude of turbulence as a function of position 
in the disk. Lee & Draine (1985) discuss the effect of 
a random magnetic field component on the strength of 
the observed polarization signature, and note that the 
strength of polarized emission will be reduced by a fac- 
tor F = 3/2 ((cos^ 0) — 1/3), where 6 is the angle between 
the local magnetic field and the direction of the ordered 
global magnetic field. This quantity varies from one (per- 
fectly ordered field; (cos^ ^) = 1) to zero (perfectly ran- 
dom field; (cos^ 6>) = 1/3), but the exact value depends 
on the details of the local magnetic field geometry. If 
magnetic field tangling were the sole factor responsible 
for the difference between the fiducial modeling predic- 
tion and the SMA upper limits, we would constrain F 
to be less than ~ 0.1 for the case of HD 163296, imply- 
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ing ^cos^ ^) < 0.4, which indicates an almost completely 
random magnetic field structure. 

It should also be noted that grain alignment efficiency 
would play a similar role, quantified in exactly the same 
way as F above, with indicating the angle between the 
long axis of the grain rather than the angle between the 
local and global magnetic fields (Greenberg 1968; Lee & 
Draine 1985). The Cho & Lazarian (2007) code assumes 
100% efficient alignment in regions that meet the grain 
alignment criterion (Section 4.2.2). In order to account 
fully for the suppression of polarized emission relative to 
the fiducial model, the alignment efficiency would have 
to be quite low, less than 10% in the case of HD 163296. 

The tentative single-dish detections appear to indicate 
a toroidal magnetic field geometry for the disks around 
DG Tau and GM Aur, consistent with observations indi- 
cating a dominant toroidal component to the magnetic 
field in the flattened structures around YSOs at ear- 
lier evolutionary stages (see Wright 2007, and references 
therein). However, it is also possible that the field could 
be poloidal: as discussed e.g. in Shu et al. (2007), a mag- 
netic field gathered from the interstellar medium that 
threads vertically through the disk might be expected 
to remain poloidal in geometry as it interacts with disk 
material. While the SMA limits are unable to constrain 
the magnetic field geometry, a poloidal geometry might 
be expected to reduce the expected polarization signa- 
ture particularly for the case of a face-on viewing geom- 
etry as in the case of the TW Hya disk. The effects of 
a poloidal geometry for a disk viewed at intermediate 
inclination, like HD 163296, are less clear and are not 
investigated in the context of the Cho & Lazarian (2007) 
models, although it is plausible that the strength of po- 
larized emission from a toroidal or poloidal field would 
be comparable. 

4.3.3. Scattering 

Cho & Lazarian (2007) argue that scattering con- 
tributes significantly less than thermal emission to the 
polarized ffux at millimeter wavelengths in the disk. In 
order to estimate the relative contribution of scattering 
and emission at a range of radii throughout the disk, they 
compare the product JA^scatt, where Jx is the mean radi- 
ation field and /^scatt is the mass scattering coefficient, to 
the product ^A^^abs, where Bx is the intensity of black- 
body radiation in the region of interest and /^abs is the 
mass absorption coefficient. They show that in the outer 
disk, where R > 10 AU, the ratio of these products falls 
below one (and ultimately below 0.5), indicating that 
emission is dominant over scattering in the outer disk. 

It is of interest, however, that pure scattering of light 
from a central source off of large grains in the outer disk 
should produce a polarization signal precisely orthogo- 
nal to that expected for elongated grains aligned with 
a toroidal magnetic field. While the radiation field at 
850 /im is dominated by the local conditions rather than 
a central source, as discussed in Section 4.2.2 there will 
be an overall radially anisotropic component of the radi- 
ation field that might be expected to produce a weaker, 
but still orthogonal on average, scattering signal. The 
contribution from scattering would therefore generally 
act to cancel the expected polarization signal from emis- 
sion. An estimate of the magnitude of the scattered light 
signal compared with the predicted strength of polarized 



emission is beyond the scope of this paper, but we note 
that for scattering to be the dominant mechanism sup- 
pressing the expected polarization signal, the intensity 
of polarized emission arising from scattering and emis- 
sion would have to be precisely equivalent, to within 10- 
15%, in both disks. Furthermore, since the scattering 
and emission have different wavelength dependences, the 
coincidental canceling of the emission signal would only 
occur at the wavelength of observation. In the absence of 
any expectation that these quantities should be related, 
this seems an unlikely coincidence. 

5. SUMMARY AND CONCLUSIONS 

Despite the expectation of a 2-3% polarization frac- 
tion in circumstellar disks based on previous observa- 
tional and theoretical work (Tamura et al. 1999; Cho & 
Lazarian 2007), the SMA polarimeter observations pre- 
sented here show no polarization from the disks around 
two nearby stars. With these observations we place a 
3(7 upper limit on the integrated polarization fraction of 
less than 1% and rule out the fiducial Cho & Lazarian 
(2007) models at the ^lOcr level. These represent the 
most stringent limits to date on the magnitude of sub- 
millimeter polarized emission from circumstellar disks. 
We are therefore left with the question of which model 
assumptions are unrealistic enough to account for an ap- 
proximately order-of-magnitude (at minimum) overpre- 
diction of the polarization signal from these disks. 

Among the model parameters and additional effects 
considered in Section 4, several seem unlikely as the 
source of the suppression of polarized emission. The 
critical magnetic field strength expected for alignment 
seems reasonable relative to theoretical expectations and 
observations. An almost completely random magnetic 
field with no dominant toroidal (or poloidal) component 
would also be surprising, although a poloidal field geom- 
etry would be expected to significantly weaken the po- 
larized emission arising from a face-on disk like TW Hya. 
Scattering is expected to be weak, but it should produce 
a polarization signature perpendicular to that expected 
for emission from aligned grains. However, scattering 
and emission signals would have to cancel nearly per- 
fectly in order to account entirely for the low observed 
polarization fraction. Nevertheless, there are promising 
candidates to describe how the suppression of polarized 
emission might have occurred. Cho & Lazarian (2007) 
assume 100% efficient alignment of grains with the mag- 
netic field in regions of the disk where the alignment cri- 
terion is met, which is overly optimistic and now known 
to be unrealistic (see discussion in Section 4.2.2). In light 
of the recent work on the quantitative theory of grain 
alignment (Lazarian & Hoang 2007; Hoang & Lazarian 
2008, 2009), the Cho & Lazarian (2007) result may be 
considered an upper limit to the theoretical expectation 
for the polarization properties of disks. A reduction to 
10% efficiency, which is within the expectations based 
on recent developments in grain alignment theory, could 
alone explain the low polarization fraction observed. An- 
other possibility is that the grains contributing most of 
the polarized emission in the model are well (or not so 
well) aligned, but rounder than the cross section ratio as- 
sumed in the initial model and therefore inefficient emit- 
ters of polarized radiation. This is also reasonable based 
on a rough estimate of the timescales and shapes ex- 
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pected for collisional growth of elongated ISM grains. 

While each of these factors would have to be substan- 
tially different from what is expected in the initial model 
to alone account for the low polarization fraction, it is 
of course entirely possible that several effects are play- 
ing a combined role. For example, grains with a cross 
section ratio of 1.5 instead of 2.1 could combine with a 
50% alignment efficiency to account entirely for the dif- 
ference between observations and models. A small degree 
of field tangling (expected because of turbulence in the 
disk) could further reduce the expected polarization sig- 
nature. While we cannot constrain precisely which fac- 
tors are contributing in which proportions to the suppres- 
sion of polarization in the disks observed with the SMA, 
we identify these three factors (grain elongation, align- 
ment efficiency, and field tangling) as the most plausible 
sources of the suppression of polarized emission. They 
produce the greatest change in polarization properties 
within a reasonable range of parameter values, and there 
exists a theoretical justification for why they should exist, 
even if the magnitude of the effect is not well constrained. 

Future observations with higher sensitivity may be able 



to disentangle these effects to some extent, particularly 
the degree of field tangling. It would also be useful to 
obtain high spatial resolution observations of the disks 
with tentative detections of a 2-3% polarization fractions 
to confirm the strength and origin of the emission on 
small spatial scales, and to expand the sample size in 
order to determine whether the low polarization fraction 
constrained by the SMA is universal for disks around 
young stars. 
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